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ABSTRACT 

The few known 7-ray binary systems are all associated with variable radio and X- 
ray emission. The TeV source HESS J0632+057, apparently associated with the Be 
star MWC 148, is plausibly a new member of this class. Following the identification 
i of a variable X-ray counterpart to the TeV source we conducted GMRT and VLA 

^T) • observations in June-September 2008 to search for the radio counterpart of this object. 

A point-like radio source at the position of the star is detected in both 1280 MHz 
GMRT and 5 GHz VLA observations, with an average spectral index, a, of ^0.6. 
In the VLA data there is significant fiux variability on '^month timescales around 
I the mean flux density of wO.3 mJy. These radio properties (and the overall spectral 

energy distribution) are consistent with an interpretation of HESS J0632-I-057 as a 
lower power analogue of the established 7-ray binary systems. 
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1 INTRODUCTION 

There are three firmly established >GeV em itting bi- 
naries; PSRB1259-63 l|Aharonian et al.l l2005bl), LS 5039 
llAharonian et all [2005al ') and LSI-l-61303 (|Albert et all 
I2OO6I ). all systems composed of a compact object and a high- 
mass star. LS 5039 harbours an 06.5V type star in a 3.9 day 
orbit, the massive star in the other systems is of Be type. 
PSRB1259— 63 is the only one of these systems in which the 
nature of the compact companion (a 48 ms period radio pul- 
sar in a 3.4 year orbit around the star) has been identified. 
All three systems display similarities in their radio to TeV 
7-ray spectral energy distributions (SEDs) including a hard 
X-ray spectrum and a softer TeV spectrum, radio emission 
above 1 GHz and variability in all bands. There is also evi- 
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dence for a si ngle TeV flare from the luminous X-ray binary 
Cygnus X-1 (| Albert et al.ll2007l ). a system of a black hole 
and an 09.7 lab star in a 5.6 day circular orbit. 

The p oint-like TeV source 

HESS J0632-I-057 jAharonian et al.l bOOTl ). has been 
suggested as a new member of this class of objects. 
This suggestion was based on associations with the mas- 
sive star MWC148 (HD 259440) and the unidentified 
ROSAT and EGRET sources IRXS J063258.3+054857 and 
3EG J0634H-0521, and would imply that MWC148 has 
a hidden compact companion. X-ray observations of this 
source with XMM-Newton in 2007 revealed a point-like, 
hard spect rum (r~1.26). X-r ay source at the position of 
MWC148 iHinton et al.ll2009h . Significant variability was 
detected on timescales of a few hours. The association of 
this X-ray source with the TeV source (and both sources 
with MWC 148) seems highly likely and the spectral 
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Figure 1. Images of the region surrounding HESS J0632+057 (shown as a black square with positional error bars) obtained with the 
GMRT at 1.28 GHz (left panel) and the VLA in D-configuration at 5 GHz (right panel). Five point-like sources are detected in both 
images (note that the increased angular resolution of the GMRT observations show that source #4 is in fact a double source). Only 
one source (source #3) lies within the rms size limit of the TeV emission from HESS J0632+057 (dashed black circle). The sources are 
labelled in order of decreasing R.A. The flux density contours are at 0.4mjy/beam and 2.0mjy/beam in the GMRT map (beamsize 2") 
and at 0.2mjy/beam and l.OmJy/beam in the VLA map (beamsize 14"). 



properties of the X-ray/7-ray source are consistent with 
the known 7-ray binaries. Very recently the VERITAS 
collaboration have published flux upper limit s which imply 
varia bility in the TeV emission of this object (|Acciari et al.l 
|2009| ). 

Historical studies of MWC 1 48 have determined the 
star to be of spectral type BOpe (iMorgan et al.lll955l ). Be 
stars are classified by the presence of Balmer emission lines. 
MWC 148 was rejected as a Herbi g AeBe star due to a non- 
detection during the IRAS survey (|The et al.lll994h and thus 
is believed to represent one of the 'classical' Be stars. These 
stars are fast rotators, though t to be spinning at 5 0-90% 
of their critical velocity (note iGutierrez-Soto et al.l (|2007l ) 
give V sm I of MWC 148 to be 430 kms"^) with a large stel- 
lar wind and a high mass loss rate. The fiattening of the 
circumstellar envelope produces a global polarisation of the 
envelope of a few percent. The most recent measured values 
of th e polarisation of MWC 148 are 3-4% (|Yudin fc Evans! 
1 19981 ). The star also exhibits a s ignificant infrared excess 
(see 2MASS. ISkrutskie et al.ll200(j ) consistent with the exis- 
tence of an extended envelope. The source lies between the 
Monoceros Loop supernova remnant and the star forming 
regions of the Rosetta Nebula (d «1.4kpc). A compatible 
distance estimate of Ril.Skpc is calculated using the appar- 
ent visual magnitude of MWC 148, Mv =9.1 l|H0g et all 



|2000|) . as described in iHinton et all l|2009l ). No variability 
was detected from the star duri ng the ASAS-3 optical sur- 
vey (|Gutierrez-Soto et al.l |2007| ). Bp stars are defined by 
their u nusual surface abundances and strong (up to kG, 
see e.g. Borra fc Landstreetl Il979l ) surface magnetic fields 



(|Phillips fc Lestrade 19881 ). The confinement of the stellar 



wind in these strong fields is thought to lead to strong shock 
heating and variable X-ray emission. iTownsend et al.l (|2007l ) 
speculate that such a process may be capable of accelerating 
particles up to TeV energies. 



All of the known binaries have been dete cted at radio 
wavelengths at fiux levels of 10-100 mjy (see lDubusll2006l . 
and references therein). The radio spectra of LS5039 and 
LSI -1-61 303 have spectral indices close to a = 0.5 (where 
oc V"). In LSI +61 303 the radio flux is modu- 
lated with the orbital period, whereas in LS 5039, 
the radio emission (unlike the X-ray and TeV emis- 
sion) is not (see e.g. Clark et al. 2001). No targeted 
radio observations of MWC 148 had previously been con- 
ducted. Here we present new radio observations resulting 
in a detection of the system with both the VLA and the 
GMRT. 
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Figure 2. XMM-Newton count map (see iHinton et al] |2009| . 

for details) overlaid with VLA 5 GHz contours (0.2 and 0.4 
mjy/beam; red) and a single GMRT 1.28 GHz contour (0.2 
mjy/beam; green). Beam sizes as in Fig. 1. The position of 
MWC 148 is shown with a white star. The HESS best-fit posi- 
tion (black square) and Itr errors are also shown. 

2 THE NEW RADIO DATA 
2.1 VLA data 

HESS J0632-I-057 was observed for 6 hours at 5 GHz with 
the Very Large ArrajQ (VLA) in D-configuration (angu- 
lar resolution Ril4") during July - September 2008 (pro- 
gram AS944). The time was divided into three 2 hour ob- 
servations separated by ~1 month. The data were cali- 
brated using the NRAO AIPS (Astron omical Image Pro- 
cessing System; iBridle fc G reiscn 1994) software pac kage 
and then loaded into DJfM^P ( Sheph erd et al.l Il994l) for 
additional editing and imaging. The flux density scale was 
set using a scan of 3C 147 at the end of each observation 
and the phase was monitored with scans of the calibration 
source 0632-1-103. The off-source rms in each observation 
(0.03 - 0.04 mJy/beam) was estimated from large source- 
less boxes, far from the phase centre and was found to be 
close to the thermal noise limit for a 2 hour observation of 
«0.03 mJy/beam. 

Five unresolved sources were found above 5(t within 
the 9' primary beam radius of the best fit position of 
HESS J0632-I-057. One of these sources (hereafter 
source #3) is located at 06''32'"59.24'' ± 0.3", 
+05°48'00.8" ± 0.3", consistent with the position 
of MWC 148. This position was derived from a two- 
dimensional Gaussian fit to the brightest single VLA 
observation. The error on this measurement is sta- 
tistical only. Emission is detected from this source in all 
three observations and the measured flux varies significantly 
from observation to observation: from 0.19 ± 0.04 mjy to 



0.41 ± 0.04 mJy (x^/dof for a constant fit = 19.4/2, chance 
probability — 6.1x10"^), see Fig. 3. Additionally, since 
each 2 hour observation consisted of four 21 minute scans 
of the region surrounding HESS J0632-I-057, we searched for 
shorter timescale (intrahour) variability. No evidence was 
found for short-term (~1 hour) variability of source #3 in 
the scan-by-scan light-curves. Source #3 was modelled with 
an elliptical 2D Gaussian in the map plane resulting in a 
size of 7" by 4" (la), consistent with an unresolved source. 

A further 3 hour observation at 5 and 8.5 GHz with 
the VLA in the high resolution A-configuration (program 
AS967) was taken in October 2008. However, source #3 was 
not detected during this observation, presumably due to a 
low flux state of this object. A plausible alternative reason 
for the non-detection is that the source is extended on scales 
signiflcantly larger than the 0.4" beam (but smaller than 
the 2" 1.28 GHz GMRT beam), reducing the signal/noise 
achievable in this configuration. To circumvent this issue, 
the A configuration data were tapered to match the 2" 
1.28 GHz GMRT beam and rms point source limits at 5 and 
8 GHz were measured and plotted in Fig. 3. All four of the 
field sources detected during the D configuration observa- 
tions were also present in the A configuration observations. 

The next closest source to MWC 148 (source 
#2) lies on the edge of the rms size limit of the TeV 
emission, at 06''33™06.01= ± 0.3", +05°47'49.63"± 
0.3". The source displays a power-law spectrum with 
Q = 0.9 ± 0.1 and has a flux density at 5 GHz of 0.7 ± 
0.2 mJy/beam. A multi-wavelength catalogue search 
resulted in no counterparts at this position. One of the 
field sources (source #1) shows significant variability both 
on the '^month separation timescale of these observations 
and on shorter (^hour) timescales. This source is located 
at 06h33m06.59s ± 0.3", -H05d53'39"2 4 ± 0.3", and is visi - 
ble both in the NVSS archive (1.4 GHz. lCondon et al.lll998D 
with a fiux of ~15 mJy (cf ~12 mJy in the GMRT 1.28 GHz 
observations) and in the XMM-Newton image of this field. 
Constant fluxes were recorded for all remaining field sources 
in all observations. 



2.2 GMRT data 

We also observed the fleld of HESS J0632-I- 057, using the 
Giant Metrewave Radio Telescops0 (GMRT) l|Swarup et al.l 
Il99ll ) during June - September 2008. The GMRT array has 
30 antennas, arranged in roughly a 'Y' conflguration over 
25 km area. The target field was observed for a total of 36 
hours with 64-6 slots of 3 hours each at 1280 and 610/235 
MHz. We observed simultaneously at 235 and 610 MHz, us- 
ing synthesised bandwidths of 6 and 16 MHz, respectively. 
The 1280 MHz observation was carried out with a band- 
width of 16 MHz in each of the two available sidebands. 
Each frequency channel is 125 kHz in width, enabling the 
removal of narrow-band radio frequency interference (RFI). 
The sources 3C 147 and 3C 48 were observed at the begin- 
ning and end of the observations and used as amplitude 
and bandpass calibrators to set the flux density scale. The 



^ The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 



^ GMRT is run by the National Centre for Radio Astrophysics 
of the Tata Institute of Fundamental Research. 
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sources 0632+103 and 0521+166 were used as phase calibra- 
tors at 1280 and 610/235 MHz respectively. All frequency 
channels affected by strong narrow-band RFI were flagged 
and removed from the entire data set using the AIPS soft- 
ware package. 

A faint, unresolved radio source is detected at the posi- 
tion of MWC 148 (see Fig. 1) in five of the sbc 1280 MHz 
observations, with an average flux of 0.68 + 0.04 mjy 
(see Fig. 3). The best fit position of this source 
from the full data set is 06''32™59.2r + 0.5", 
+05°48'01.2" ±0.5". There is no significant evidence for 
variability in these data over the week to month timescales 
of the observations. MWC 148 was not detected above the 
3o" level in the lower frequency observations. Upper limits 
on the flux at 610 MHz and 235 MHz have been calculated 
and included on the SED shown in Fig. 4. The typical noise 
levels at 610 and 235 MHz were 0.74 and 1.82 mjy/beam 
respectively, the <lGHz data were more strongly affected 
by RFI, resulting in a substantial increase in the rms noise. 



3 DISCUSSION 

The radio properties of source #3 can be summarised 
as follows. The best fit positions of source #3 measured 
with the VLA and GMRT are consistent both with each 
other and with the position of MWC 148 (06*' 32^59.236", 
+O5°48'Ol.O4": lH0g et alll2000l ). The new GMRT data pro- 
vides the best constraint on the size of the radio emitting re- 
gion in MWC 148. The upper limit on the source size (along 
the short axis of the GMRT beam) is 2". Assuming a dis- 
tance to the star of 1.5 kpc, this corresponds to an upper 
limit on the emission region of 3000 AU. Significant variabil- 
ity was detected in the 5 GHz emission from MWC 148 at 
the 5cr level. Variability at the same level may be present in 
the 1.28 GHz data but not detectable above the noise. The 
non-detection during the VLA A configuration observations 
further demonstrates the variable nature of the source. The 
variability timescale must be longer than the Ri2hour scans 
and shorter than the ~month timescale separation of the 
observations. An estimate of the spectral index, based on 
the average fiux at 1.28 GHz and 5 GHz, was calculated for 
source #3 to be ar = 0.6 ± 0.2. However, this value should 
be treated with caution due to the non-simultaneous nature 
of the 1.28 GHz and 5 GHz observations. 

The new radio source l ies less than ^1" f rom the cen- 
troid of the X-ray emission (iHinton et al.|[2009l ). and within 
the la error box of HESS J0632+057 (see Fig. 2). Thus 
we confidently identify this emission as the radio counter- 
part to MWC 148. It seems very likely that all these ob- 
jects are associated and therefore we can create an SED as 
shown in Fi g. 4. The re cent publication of the Fermi Bright 
Source List (lAbdo et al . 2009) allows us to place constraints 
on the GeV emission from HESS J0632+057. As no source 
was detected (above lOo" in the first 3 months of operation) 
at the position of HESS J0632+057, the EGRET source 
3EG J0634+0521 can be ruled out as a potential counter- 
part. Limits for the GeV emission after 3 months, and 1 
year of Fermi data have been added to the SED. 

The detection of TeV emission implies particle acceler- 
ation is taking place in the source. The most natural ex- 
planation for the observed radio emission is then optically 



thin synchrotron emission from accelerated electrons. Syn- 
chrotron radiation with a spectral index ~0.5 implies an 
underlying electron spectrum with dN/dE oc . This can 
be interpreted as the injection spectrum. IC cooling in the 
Thomson regime is hard to avoid for these electrons. If the 
observed spectrum is cooled, then this implies a very hard 
injection spectrum. This scenario was the one presented in 
iHinton et all l|2009l '). in which there is a low energy cut-off 
in the injection spectrum, resulting in an effectively mono- 
energetic injection for electrons below this cut-off energy. 
Synchrotron cooling of such an injection will produce a time 
averaged electron spectral index of 2 and an emission spec- 
trum with q=0.5. T he simple one-zone model presented in 
IHinton et all (|2009l ') and adjusted to fit the X-ray and TeV 
data, results in a radio spectrum remarkably similar to that 
observed. A better fit to the radio data can be achieved 
by adjusting the low-energy cut-off or a change in injec- 
tion index from 2.0 to 1.9 (see Fig. 4). Ass uming a one- 
zone m odel with the parameters described in IHinton et al.l 
(|2009l '). B = 70 mC, (7rad ~ lergcm"^, the radio emission 
at 5 GHz would be produced by ~90 MeV electrons having 
a characteristic cooling time due to IC losses of ~2 days. 
Similarly, 1.3 GHz emission would be produced by 50 MeV 
electrons with a cooling time of ~4 days. The observed vari- 
ability timescale, in common with that seen in X-rays, is 
comparable with the expected cooling times, even if the ra- 
dio emission region is somewhat larger than the few AU 
assumed. We note that if the emission region was consider- 
ably smaller than this both significant free-free absorption 
at 1 GHz in the stellar wind and 7-7 absorption at a few hun- 
dred GeV on stellar photons would likely become apparent 
at some orbital phases. At this stage neither effect can be 
excluded. Strictly simultaneous multi-frequency radio data 
would be required to check for radio absorption features. 

There are two main scenarios that could plausibly ex- 
plain the observed emission from HESS J0632+057; firstly 
that the star is isolated, with a magnetically confined wind 
region and secondly that the star is in a binary system with 
an unseen compact companion. Although there are some in- 
teresting similarities between MWC 148 and the archetypal 
isolated magnetic Bp star a Ori E, including similar X-ray 
variability timescales (|Skinner et al.ll2008l : lAharonian et al.l 
l2007fl . there are also significant differences. The X-ray emis- 
sion from a Ori E can be well fit by a two-temperature apec 
model with a hot component of ~2.6keV whereas the X-ray 
emission from MWC 148 can only be fit with such a model if 
the temperature of the hot component is well above 10 keV, 
due to the substantial fiux in the 5-10 keV band. At radio 
wavelengths a Ori E displays a fiat spectrum [a ~ 0)between 
2-20 GHz, whereas the emission from MWC 148 appears to 
have a negative spectral index. The radio light-curve of o Ori 
E is strongly modulated with the orbital period. This may 
also be the case with MWC 148, but further observations are 
clearly needed to characterise its radio variability. 

As can be clearly seen in Tab. 1 there is a 
strong resemblance between the known 7-ray binaries and 
HESS J0632+057. Although HESS J0632+057 is a much 
weaker source in all wavelengths, the similar spectral prop- 
erties make a convincing argument that this object does 
indeed represent a new 7-ray emitting binary system. 
The apparent variability in the TeV emission of 
HESS J0632+057, as inferred by the VERITAS col- 
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Figure 3. The radio flux density of MWC 148 as a function of time (note that JD 2454620 is the 2^1^ June 2008). The black star-shaped 
points denote observations taken with the VLA and the red circular points are from the GMRT. The dotted lines represent the weighted 
mean flux at each frequency. Flux density upper limits (at 99.8% confidence level) from the VLA A configuration observations are plotted 
in blue. 
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Table 1. Properties of the TeV emitting binaries, adapted from lPaxedesI l|2008l V Note that the spectral indices are defined by Fi, oc v~'^ 
or equivalently dN/dE oc All luminosities are in units of lO^^ergs"^ except the radio luminosities which are in units of 

l O'^^ergs"^. Lumino sit ies are given for t he fol low ing ranges: L,-: 0. 1 -100 G Hz, Lx: 1 - 10 ke V, Lgbv: 1 ~ lOGeV, Lxev: 0.2 - 10 TeV . 
^Marti et all lll998l). ^iTakahashi et al.l (|2009|), MAharonian et al.l ll2006tl. iGregorvl ll2002h -I Albert et al] 1I2OO6I) , lAlbert et ahl ll2009t). 
^ Johnston et al.l ||2005D. jEsposito et al.l (120071) . ^ lAharonian et al.l ll2005bh . ^Pandev et al.l | |2007^ . jMiller et al.l ||20()5^. ^ Hinton et al.l 
1I2OO9I '). ^ Aharonian et al 1 ll2007t) . * GeV luminosity measurements are from the Fermi Bright Source List llAbdo et al.ll2009ir with upper 
limits estimated from non-detection at a lOo" level after three months of observations. Note that the association of GeV emission with 
LS 5039 is based only on positional coincidence. No GeV detection reported yet. 



laboration (|Acciari et al.l l2009h . is also consistent 
with the properties of the known 7-ray binary sys- 
tems. All objects in Tab. 1 exhibit a hard [ax < 1) spec- 
trum in the X-ray domain, with HESS J0632+057 having a 
particularly hard spectrum. The TeV spectral indices of all 
objects are consistent with each other and all considerably 
softer than the corresponding X-ray spectrum. All objects 
appear to have their peak energy output in the MeV to GeV 
range. Note that GeV emission has only been de tected using 
Fermi coincident with LS I -f 61 303 and LS 5039 (|Abdo et al.l 
I2OO9I ). Upper limits for the other sources have been esti- 
mated using detected sources close to the objects of interest 
to determine the local detection threshold. The 7-ray bina- 
ries all display variable radio emission with a wide range of 
spectral properties including lo w energy spectral turnovers 
(see e.g. iGodambe et al.l l2008h . The clear outlier in this 
group is Cygnus X-1, for example, although the radio and 
VHE 7-ray luminosities are similar, the X-ray luminosity of 
Cygnus X-1 is orders of magnitude greater than the other 
objects. Resolved radio emission has been detected in both 



LSI + 61 303 and LS 5039 (|Paredes et al.l Il998l : iRibo et all 
|2008| ) with emission regions of the order of ~10AU. Our 
current observations are unable to probe such small scales 
and more data is required in order to resolve or further con- 
strain the size of the emitting region and move towards a 
better understanding of this source. Whilst the detection of 
variable radio emission from MWC 148 further strengthens 
the case that HESS J0632+057 is indeed a TeV binary, the 
scenario that this object is an (extremely unusual) isolated 
massive star remains an intriguing possibility. 
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